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ABSTRACT 


The  GalerkiD  method  is  applied  to  Flugge’s  differential  equations  for  the  vibration  of  a 
cylindrical  shell  under  axial  load  and  external  pressure  to  obtain  a  3N  x  3N  characteristic 
equation  in  matrix  form.  N  is  the  number  of  terms  in  the  assumed  series  of  displacement 
functions  for  the  uf  v,  and  w  displacements  which  can  be  selected  to  satisfy  various  boundary 
conditions.  For  the  freely-supported  cylinder,  an  exact  solution  exists,  and  the  various 
assumed  modes  uncouple  reducing  the  problem  to  the  solution  of  a  3  x  3  characteristic 
equation  for  each  mode. 

The  third  order  characteristic  equation  for  the  freely-supported  cylinder  was  solved  for  a 
wide  range  of  shell  parameters.  The  natural  vibration  frequencies  and  buckling  values  for 
axial  load  and  external  pressure  for  all  three  eigenvalues  associated  with  each  mode  are 
presented  in  a  series  of  figures.  The  square  of  the  vibration  frequency  for  any  mode  was 
found  to  vary  linearly  with  axial  load,  and  approximately  linearly  with  an  external  pressure 
loading  for  modes  with  two  or  more  circumferential  waves. 
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SECTION  I 
INTRODUCTION 


In  recent  years  thin  shells  have  been  used  extensively  for  aerospace  applications  because 
of  the  high  structural  efficiency  of  this  type  of  construction.  Space  capsules,  boosters,  space 
stations,  and  large  fuel  tanks  for  hypersonic  cruise  vehicles  are  examples  of  stiffened  and 
unstiffened  shell  configurations.  In  many  cases  these  shells  are  pressurized  and  subjected 
to  some  type  of  loading.  These  loads  and  internal  pressure  can  have  a  significant  effect  on 
the  vibration  characteristics.  The  purpose  of  tnis  investigation  is  to  examine  the  combined 
effects  of  axial  load  and  external  pressurization  on  the  *  ibration  characteristics  of  thin 
cylindrical  shells  and  present  the  results  in  a  suitable  form  for  use  in  design. 

Arnold  and  Warburton  (References  1  and  2)  examined  the  vibration  characteristics  of 
freely-supported  shells  and  plotted  results  for  several  radius  to  thickness  ratios.  They  in¬ 
dicated  that  three  natural  frequencies  exist  for  each  modal  pattern.  The  different  frequencies 
are  associated  with  motions  that  are  primarily  radial,  longitudinal,  or  torsional.  The  lowest 
frequency  is  usually  associated  with  motion  which  is  primarily  radial.  Forsberg  (Reference  5) 
obtained  an  exact  solution  for  the  cylindrical  shell  using  Flugge’s  equations  of  motion  with 
various  homogeneous  boundary  conditions.  Fung  (References  7,8,  and  10)  obtained  analytical 
and  experimental  results  for  the  effects  of  internal  pressurization  on  the  vibration  character¬ 
istics  of  cylindrical  shells.  Fung  and  Sechler  (Reference  9)  and  Weingarten,  Morgan,  and 
Seide  (Reference  13)  obtained  results  for  the  buckling  of  cylindrical  shells  subjected  to  in¬ 
ternal  pressurization  and  axial  load.  A  recent  paper  by  Herrmann  and  Shaw  (Reference  11) 
on  the  vibration  of  thin  shells  under  initial  stress  included  the  change  in  magnitude  and 
direction  of  the  applied  load  as  a  result  of  deformation  and  contains  a  comparison  of 
experimental  data  witn  various  shell  theories. 

This  investigation  applies  the  Galerkin  method  to  Flugge’s  cylindrical  shell  equations  to 
formulate  a  characteristic  equation  valid  for  arbitrary  boundary  conditions.  By  selecting  a 
series  of  functions  for  jach  of  the  displacements  u,  v,  and  w  various  boundary  conditions  caa 
be  handled.  For  the  special  case  of  an  unloaded,  freely-supported  cylinder,  Arnold  and 
Warburton  (Reference  1)  obtained  an  exact  solution.  This  same  solution  was  obtained  by 
Flugge  (Reference-  4)  for  the  buckling  of  a  cylinder  subjected  to  axial  load  and  internal  pres¬ 
surization,  and  used  by  Herrmann  (Reference  11)  indicating  that  the  same  deflection  mode 
shape  satisfies  the  combined  problem.  When  a  series  of  functions  of  this  form  were  substituted 
into  the  equations  obtained  in  this  investigation,  it  is  shown  that  the  terms  in  the  series  un¬ 
couple  and  Flugge’s  buckling  results  are  obtained  if  the  frequency  is  assumed  to  be  zero. 
Vibration  results  comparable  to  those  of  Arnold  and  Warburton  are  obtained  if  the  load  pa¬ 
rameters  are  assumed  to  be  zero. 

For  this  special  case  the  axial  load  parameter  and  frequency  parameter  can  be  combined 
into  a  single  eigenvalue  expression.  The  eigenvalue  spectrum  for  this  case  is  plotted  in  non- 
dimensional  form  for  a  wide  range  of  shell  parameters  covering  radial,  longitudinal,  and 
axial  deformation  modes.  From  these  charts  and  tables  the  critical  buckling  conditions  or 
vibration  frequencies  can  be  obtained  for  the  freely-supported  cylinder  subjected  to  an  axial 
load  and  external  pressurization. 
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SECTION  II 

THEORETICAL  DEVELOPMENT 


The  coordinate  system  used  in  the  analysis  is  the  same  as  that  of  Flugge  and  is  given  in 
Figure  1.  The  equations  of  motion  for  thin  cylindrical  shells  including  the  effects  of  a  uniform 
normal  pressure  and  an  axial  load  can  he  obtained  by  adding  the  inertia  terms  to  Flugge’s 
static  equations  given  in  Reference  4. 
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The  equations  of  motion  can  be  written  in  the  following  matrix  form 
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The  differential  operators  L. .  can  be  written  as  follows  after  multiplying  Equations  1  and  2 
by  -1. 
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dr 


(8) 

(9) 

(10) 


.  .  t  i  .  »  d  ,  k  ( 3 -v)  2  d 

Lz3l  _{1_q!  >_  4-— -  a 


d<p 


d>*d<t> 


HD 


31 


,  .  x  d  .  k  (l-l/).  od 

a  (j/  4  q,  )—  4  — — 


t32=  11  q'  }  djy 


dxd<p 

k(  3  - v  )  oV 
T 


k  o  3  d3 
dx3 


r  cV 

dx4  dx2d<^>~ 


LJ3:  *  ■**  k 


2a  2  d4 


dx  dcf> 

d4 


4-M_  4,  I 
d<£4  d<£2  J 


(12) 


+  0, 


d<£‘ 


+-  q,  a 


2  d 


at 


(.4) 


To  obtain  a  solution  by  Galerkin’s  method,  the  equations  of  motion  are  multiplied  by  the 
variation  of  the  displacements^  x  ,  and  the  result  integrated  over  the  volume  r  f  the  shell  (Ref¬ 
erence  12).  Since  the  integration  over  the  thickness  has  already  been  performed  in  deriving 
the  equations,  none  of  the  factors  is  a  function  of  r .  It  is  therefore  sufficient  to  integrate  over 
the  shell  area.  This  can  be  written  as  follows: 


IJ  {Sx}  [  l]  {  x }  «  =o 


( i5) 


The  displacements  u,  v,  and  w  are  periodic  in  cf>  with  period  2v .  Simple  harmonic  motion  is 
assumed  and  the  displacements  are  written  in  the  form  of  a  product. 


u  =  U  (  x  )  co  s 


(  16) 
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Substituting  Su»  8v,  and  8w  from  Equation  23  into  Equation  24, 

//Jcosn^ewt  I  SAj  <t>j  (x  )  [  L„  u  fL|2  v  F  L,3w] 


Hx)  |  ^  y  r 

+  s.nn(/>e  £  Bb.$.  <  *  >  [  L  2  ,  u  +  L  2  2  v  +  L  23  *  1 

i  =  I  •* 

N 

+  cosn<£e'W  ZSC,  C  .  >  [Lj,u+  LJ2  »  f  L3J  w  ]  ) 


dA  =0 


Since  8A.,  SB.,  and  8C.  are  arbitrary,  the  only  way  that  the  above  equation  can  be  identically 


zero  is  that  each  integral  vanish  individually  (Reference  12), 

J  J  cosn(/>el<iJl  SAj  (x)  [  L(|  4j  +  L|2v  +  L|3w  ldA-0 

L  J 


d  A=  0 

(26) 

i  ;  |  O 
'  .  *-  » 

•  N 

dA  --  0 

(27) 

'  =1,2, 

N 

dA  =0 

(28) 

i  =  1,2, 

N 

JJ  cosn^eiajt  8C(  (  x)  [  L  ,,  u  + LJ2  v  +  LJ3  w  ]  dA^O  (28) 

1=1,2,  N 

Since  SA.,  83.,  8C.f  and  elwt  are  independent  of  the  coordinates  of  dA,  these  can  be  fatten 
out  of  the  above  integrals,  to  obtain 

J J  cos  n  c£>  ( x )  [L(|  ii-Fl-!  v  +-  L  |3  w  j  dA  =  0  (29) 

A  '  ”  *‘  =  I,  2  ,  •  •  N 

JJ  si  n  n  <^>  (*)  |i  2(  o  f  L22v  f  L  2j  *  ]  dA  =  0  (30) 

A  '=1,2,  N 

JJ  cos  n  (  X  )  [l  3,  u4-L3Zv+L33  w  ]dA=0  (31) 


t2|u-tL22v4-L23w  j  dA  =  0 


.=  1,2,  N 


i  =  i,  2,  -  N 


The  expressions  for  u,  v,  and  w  given  by  Equations  16  to  18  are  substituted  into  Equations  29 
to  31  with  dA  =  ad<£  dx  and  the  integration  performed  on  From  the  orthogonality  of  sin 
n<£  and  cos  n<p  over  the  interval  0  to  27r, 


t—  M 

J  sinZn0d^)  =  7r 


2tt  z 

J  cos  n  <^>  d  -  it 
0 


o 
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U,  V,  and  W  as  given  by  Equations  19  to  21  are  substituted  into  Equations  38  to  46.  These 
expressions  are  substituted  into  Equations  35  to  37  with  N  =  M  to  obtain  the  following  form 
for  the  characteristic  equation. 


(47) 


where  the  submatrices  K1^  and  M1'*  can  be  written 


0 

0 

K,J  = 

< 

o;1. 

0 

M22 

0 

-  °! 

0 

0 

m" 

M33 

_ 

The  terms  in  the  D^,  and  Y.  submatrices  are. 


D'J  =  f  <t>"  (x)  [ 

12  J  I  L 


l 

dx* 

a  (!+  J/)  d  1 


n  -  1  $  (x  )  dx 
dx  J  i 


°,3  ''I  U,di 

0  dx 


ii  f  * .  v  r  a(!  +^)n  d  u 

°2,  -I  U)L  Y~  17  J  *> 


dx 


P  2 

jj  r  .v  r2  z  r(i-i') 

°ZZZ5  %  <x)  C-q.J-a  (l  +  3k) -qg  ]—  }  $  U  )d« 


-  o2( 3- «/ )  d 


Dtt*/ 


d  1  * 

kn  — =-  <P.  (  x  )  dx 

2  dx2  J  J 


Dd,  * l  **'{■['♦%  »]~-o5k  !,)  a. 


ij  rfw  r 
D*z  i (x) 


o  n(3 -i/  )  d 

- k  — - 

2  dx* 


i*; 


(  x  )  dx 


(48) 


(49) 

(50) 

(51 ) 

(52) 

(53) 

(54) 

(55) 

(56) 
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A  solution  is  obtained  by  selecting  a  series  of  displacement  functions  $>u(x),  (£>'  (x),  and 

$w(x)  for  U(x),  V{x),  and  W(x)  respectively  which  satisfy  the  appropriate  boundary  conditions 

for  the  shell.  Integrals  (Equations  49  to  61)  are  evaluated  a"d  the  various  K1^  and  M1J  matrices 
given  by  Equation  48  are  formed  and  substituted  into  Equation  47  which  can  bo  solved  by  stan¬ 
dard  eigenvalue  techniques. 

If  the  series  of  displacement  functions  selected  for  $>Uvx),  3>v(x),  and  'PW(x)  each  form  an 
orthogonal  set,  the  generalized  mass  matrix  M1-*  in  Equation  47  will  be  diagonal. 

One  set  of  functions  which  can  be  used  for  the  series  of  assumed  displacements  are  the 
beam  functions,  whicn  satisfy  the  following  differential  equation: 

d 

- —  -  k  <£  -  0  (62) 

d  x4 

The  functions,  <$> ,  taken  as  the  assumed  modes,  satisfy  various  combinations  of  the  following 
boundary  conditions  at  x  =  0,  and  x  =  £ 


<t>  : 

0 

(63) 

0 

1 64) 

<E>"  = 

0 

(65) 

0 

(66  J 

The  beam  functions  satisfying  Equation  62  for  various  combinations  of  the  boundary  con¬ 
ditions  (Equations  63  to  66)  are  orthogonal.  Integrals  of  the  beam  functions  and  their  deriv¬ 
atives  have  been  evaluated  by  Felgar,  and  formulas  for  these  integrals  are  given  in  Ref¬ 
erence  3.  In  general,  the  terms  in  the  matrix  Klj  given  by  Equations  49  to  57  involve  integrals 

of  combinations  of  and  which  may  tie  different  functions.  These  are  not 

tabulated.  One  approach  would  tie  to  evaluate  the  various  integrals  in  Equations  49  to  57  numer¬ 
ically.  '’’’his  would  allow  an  independent  choice  of  the  function  to  represent  each  of  the  functions 
U(x),  V(x),  and  \V(x). 
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SECTION  III 

SOLUTION  FOR  THE  FREELY-SUPPORTED  CYLINDER 


An  exact  solution  can  be  obtained  for  the  vibration  of  a  cylindrical  shell  with  axial  loads 
and  internal  pressurization  for  the  following  boundary  conditions: 


w  =  0 

at 

x  =  0,  / 

!  67) 

v  =  0 

of 

X  =  0,  / 

(68) 

Mx  s  0 

at 

X  =  0,  4 

(69) 

Nx^  0 

a  t 

x  --  o ,  i 

(70! 

The  boundary  conditions  corresponding  to  Equations  67  to  70  are  satisfied  if  the  functions 
U,  V,  and  W  satisfy  the  following  conditions  at  x  -  0  and  x 


dll 


=  0 


W  - 


d  x 
V  =  0 
d2  W 


(7M 

(72) 

(73) 


dx 


These  boundary  conditions  are  satisfied  by  the  following  normalized  displacement  functions: 


;  x )  =  Jz 

C  u)  *  ST 


<Pm  Ui  =  JT 


cos 


sm 


sin 


m  t r  x 

e 

m  7 rx 

£ 

m  7T  x 

£ 


m) 

(75) 

(76) 


Tae  previous  sine  and  cosine  functions  are  orthonormal  on  the  interval  0  to  27r  and  have  the 
:o Rowing  properties: 

.£ 


/  (yr  (ST 

0  x 

f  (yrs.,i-p)(,/T 
0  c 


/■*  )(  yr 

n  - 


COS 


sin 


cos 


\ir  *  \ 

A  v 

- 1  y’ 

i  -  i 

a  x 

'  t  0, 

■  *J 

)  7 r  x  ' 

)dx 

‘  z  1 

£ 

"10, 

t*  i 

1  IT  X 

) dx 

n 

'  £ 

u 

(  77) 


(78) 


(79) 


For  this  case  the  various  integrals  in  Equations  49  to  60  are  zero  for  i  4  j  and  reduce  to  the 
following  for  i  -  j  =  m. 


mm 

D„  * 

a2  (1  -t<2) 

/  mir 

\  J 

mm 

n 

0+  v) 

l  m  ir 

0  2 

'  t 

(i-i) 


(  i  +  k  ' 


p2  £ 


(60) 


(81 ) 
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Equation  93  is  of  a  diagonal  form  and  consists  of  the  direct  sum  of  3  \  3  subnutrices.  fins 
indicates  that  an  exact  solution  exists  for  the  assumed  functions  selected.  For  a  non-trmal 
solution. 


2  2  ..H. 

( K  -  y  oj  M  ) 


,  ..22  2  2  ,.22  . 

( K  -  y  uj  M  ) 


, ^NN  2  j  NN 
(  K  -  y  to  M 


-0  (94). 


Applying  the  Laplace  expansion  to  the  above. 


or 


„ll  2  2  ,.ll  ,,22  2  2  22  J 

K  -  y  oi  M  K  —  y  uj  M 


mm  2  ?  mm 

K  -  y  ui  M  ;  0 


T  mm  ?  ?  mml  f  1 

Ik  -x  «-  m  J\ vm}  , 


..NN  2  2  ,,NK‘ 

K  -  V  u»  M 


m  :  I,  2, 


0  m  =  i ,  2 , 


=  0  '(95) 

{ 96; 
(97) 


The  expressions  for  Ku  and  M11  are  substituted  into  Equation  97.  The  resulting  expression 
is  divided  by  H,  and  y2w2  and  q9  are  taken  to  the  right-hand  side.  A  non-dimensional 

wavelength  parameter  X  is  introduced,  and  the  following  result  is  obtained: 


P  m 

m 

m 

D„ 

D 12 

0,3 

Dm 

Dm 

Qm 

21 

22 

23 

„m 

m 

L  D 

D 

D 

3i 

32 

33 

where 


,22  .2  . 
- ( y  w  - q„  X  ) 
'  2 


08) 


m  v  o 


Di i  :  x*  f  i!  +  k)  "  o  f  j  "2 


°I2'  ^21 


X  (  I  + 1/  ) 


D_  =  D 
re  3i 


u  +  q .  - 


(  !  —  v  )  2 

- ~ -  in 


-k  X3 


D22:  r>?(l  -q.  )  +  A2  f- 


'2  r  ""-‘U  ,l  +  5, 

2  J 


)] 


(  995 

( .001 

( !  0 1 ) 

<:02) 

(103) 
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D23  :U3  2  ;  (•-  q(  )  n  +  — — y — ~  kn  1 104) 

°33  =  ['  —  n2q,  -+-  k  (n2— i)£J-|-2n2KX2  +  k  A  ( 105) 

In  Equation  98  the  vibration  frequencies  and  axial  load  parameter  are  combined  to  form  a 
single  eigenvalue  expression.  This  is  discussed  in  more  detail  in  Section  IV. 

If  cu  is  assumed  to  be  zero  in  Equation  98  these  equations  are  identical  to  those  given  by 
Flugge  in  Reference  4  for  the  buckling  of  a  cylindrical  shell  subjected  to  an  axial  load  and 
internal  pressure.  If  the  axial  load  parameter  q0  and  internal  pressure  parameter  q^  are 

assumed  to  lie  zero  >n  Equ  iliun  98  these  equations  are  comparable  to  those  obtained  by  Arnold 
and  Warburton  (Reference  1).  Small  differences  which  occur  are  in  terms  containing  the 
thickness  to  radius  parameter  k  whicn is srnaUfor  cylinders  for  which  the  thin  shell  equations 
apply. 

Once  the  mode  rumbers  are  selected,  the  values  of  the  D..  terms  can  br  e/aluated  for  a  given 
shell  having  *he  assumed  boundary  conditions  and  a  given  lead  of  q^.  This  is  accomplished  by 
assigning  integer  values  to  n  and  to  m  (in  A  =  fT77--)  which  represent  the  number  of  circum¬ 
ferential  waves  and  the  axial  half-waves  respectively. 

A  solution  can  be  obtained  by  solving  Equation  98  by  matrix  iteration  or  by  reducing  Equa¬ 
tion  98  to  a  determinant  form  and  solving  the  resulting  third  order  characteristic  equation. 

j  D  —  AI  |  =0  ,  A:(/ZouZ+  q  g  £  )  (|061 

or 


Three  values  of  A  can  be  obtained  from  Equation  108  indicating  that  three  frequencies  exist 
for  each  of  the  assumed  mode  shapes.  This  result  has  been  previously  reported  by  Arnold  and 
Warburton  (Reference  2.)  and  others.  These  various  frequencies  are  associated  with  different 

A  p 

values  of  the  amplitude  ratios  ^  and  ror  one  value  of  the  frequency,  usually  the  lowest, 

C  >  A,  B  and  the  motion  is  radial  or  in  the  z  direction.  The  other  two  values  of  frequency  are 
usually  associated  with  motion  that  is  tangential,  B  >  A,  C  where  the  primary  displacement  is 
v, or  longitudinal  vibration,  where  A>  B,  C  and  the  u  displacement  predominates.  It  is  noted  that 
ior  the  cylindrical  shell  this  coupling  between  displacements  exists. 
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For  the  non-axisymmetric  modes  (n  4  o)  the  values  of  A  can  be  determined  from  Equa¬ 
tion  108  and  the  amplitude  ratios  ^  and  ~  can  be  obtained  from 


n-A 

0.2 

0  13 

12 

0Z2"A 

P23 

* 

D 

< 

1 

Q 

13 

23 

33 

D 
D 

Solving  Equation  109  for  the  amplitude  ratios  yields 


j:{ 

1  c  \ 


-o 


(109) 


_A_ 

C 


M3 

5I2 


(D22-A  )  ~D 


23 


(  D 1 1  -A)(D22-  A) 


D,Z- 


M2 


i  MO) 


"c 


Pl  3  °23 
Die 


—  ( O3  3  ~  A  ) 


0,3 

D„  -  ~  (  D  2  2  ~  A  ) 


'23 


'12 


(III) 


A  coupling  between  the  three  displacements  u,  v,  and  w  exists  for  the  non-axisymmetric  modes. 


For  the  special  case  of  axioymmetric  motion  with  n  =  o.  the  terms  D12,  D^,  D23>  and  Dj2 

in  Equation  98  are  zero,  and  the  tangential  displacement  v  can  be  uncoupled  from  the  axial 
and  radial  displacements.  With  n  =  o.  Equation  98  can  be  written  as  two  equations 


- 

D„ 


L  D 


m 

D. 


3 

m 

33 


(  c„  ) 


:(yZJ-  +Q  £  ) 


m  =  I,  2, 


1 

N 


(M2) 


[  °2  2  “  {yZ  ™  +  q2  ^  J1  0  (ll3) 

m  -  ! ,  2  , 

Setting  the  determinants  of  Equations  112  and  113  equal  to  zero  and  expanding  gives: 


A2  -  (0.,  +  D33)  A  +  D„  D33  -  D*  =0 


A  -  D 


22 


(114) 
(  115) 


Equation  115  represents  the  eigenvalue  for  pure  torsional  vibration,  and  Equation  114  cou¬ 
pled  radial  and  axial  vibration  eigenvalues.  The  amplitude  ratios  for  the  latter  are  given  by 


A 


—  D 


13 


( D , ,  -  A  ) 


(116) 
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SECTION  IV 

VIBRATION  AND  BUCKLING  RESULTS  FOR  AXIAL  LOADING 


It  was  shown  in  Section  III  that  for  the  specific  case  of  a  freely-supported  cylinder,  the 
various  terms  in  the  series  of  assumed  modes  uncouple,  indicating  that  an  oxact  solution  was 

2  2 

obtained.  The  frequency  parameter  y  u >  was  combined  with  the  product  of  the  axial  load  pa¬ 
rameter  q0  and  the  square  of  the  axial  wave  parameter  X  to  form  a  s  ingle  eigenvalue 

expression  A(X,  n). 

A  (  X  ,  n  )  -  yZ  J-  -4-  q  X2  (l!7) 


The  expressions  which  make  up  the  coefficients  of  the  characteristic  equation  contain 

the  shell  parameters  v  and  k  ,  the  wave  parameters  n  and  X,and  external  pressure  parameter 
q  .  The  parameters  u>  and  q^  appear  only  in  the  parameter  A.  When  integer  values  are  as¬ 
signed  to  the  number  of  circumferential  waves  n  and  axial  half-waves  m,  the  mode  shape  is, 
in  effect,  specified  and  solution  of  the  appropriate  characteristic  equations  yields  the 
eigenvalues  for  that  mode.  Some  typical  mode  shapes  are  shown  in  Figure  2.  It  is  noted  that 
for  the  boundary  conditions  considered,  the  eigenvalues  A(X,  n)  can  be  obtained  without  spec¬ 
ifying  a  value  for  the  axial  load  parameter  The  eigenvalues  obtained  depend  on  the  internal 

pressure,  however,  since  the  coefficients  of  the  characteristic  equations  contain  the  parameter 

V 

In  this  Section,  the  effects  of  internal  pressure  are  neglected  and  q.  is  considered  to  be 


zero.  If  the  parameter  q  is  zero,  the  eigenvalues  obtained  for  the  specified  mode  cor- 

respond  to  the  values  of  the  vibration  frequency  parameter  yV  for  the  shell  with  no  axial 
loading.  The  circular  frequency  for  the  unloaded  shell  is  ojq. 


A  (  X  ,n  )  :  y 2  (118) 

•  o 

From  the  three  values  of  A  CX,  n)  associated  with  each  mode,  the  vib:  ation  frequencies. 


a»o,  for  radial,  longitudinal,  and  torsional  vibration  having  the  specified  mode  shape  can  be 

determined  for  a  zero  value  of  axial  load.  These  results  correspond  to  those  obtained  in 
Reference  1. 


The  characteristic  equation  (Equation  108  or  Equations  114  and  115)  for  the  axisymmetric 
modes  (n  =  0)  was  solved  for  a  range  of  radius  to  thickness  ratios  from  20  to  5000.  For  large 

values  of  the  radius  to  thickness  ratio  (^-=  2000,  5000)  and  large  values  of one  of  the  roots 

was  quite  small  (A«l).  In  this  case,  greater  accuracy  was  obtained  by  solving  the  smallest 
eigenvalue  using  only  the  linear  terms  of  Equation  108.  The  results  are  given  in  Figures  3 
through  x6  for  the  first  eigenvalue  which  is  radial  deformation  (C>A,  B).  Figure  17  gives  the 
results  for  the  second  eigenvalue,  and  Figure  18  the  results  for  the  third  eigenvalue.  The 
results  for  the  axisymmetric  case  (n  =  o)  are  given  in  Figure  19. 


A  single  figure  covers  the  range  of  radius  to  thickness  ratios  from  20  to  5000  for  modes 
corresponding  to  the  second  and  third  eigenvalues.  Since  the  results  are  independent  of  the 
radius  to  thickness  ratio,  bending  effects  are  relatively  unimportant  as  compared  to  stretching 
or  membrane  effects.  The  changes  in  the  eigenvalue  spectra  with  various  values  of  radius  to 
thickness  for  radial  deformation  indicate  that  these  modes  involve  primarily  bending  effects. 
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Two  curves  are  given  for  each  radius  to  thickness  ratio  for  radial  deformation.  In  one  case 

J~K  is  plotted  versus  the  a.ial  wavelength  parameter  for  different  values  of  n  or  number 

of  circumferential  waves,  and  in  the  other  case  •/&  is  plotted  versus  n  for  different  values 
0  0 

of  — .  For  a  given  cylinder  the  ratio  —  is  fixed.  By  assigning  integer  values  to  m  in  the  pa- 
ma  a 

rameter  — —  the  values  of  for  various  modes  can  be  determined  from  the  figures..  Since 
ma  b 

is  equal  to  yw  for  q„  =  0,  the  natural  vibration  frequencies  for  the  unloaded  shell  can 
o  ^ 

be  obtained  from  these  figures.; 

For  modes  having  one  circumferential  wave  ( n  =  1)  the  amplitude  of  the  radial  and  tangential 
displacements  corresponding  to  the  lowest  eigenvalue  ore  approximately  equal  and  greater  than 

the  axial  displacement  for  values  of  ^  >  3.5  and  the  deflection  shape  is  similar  to  beam  bend¬ 
ing  with  little  deviation  from  a  circular  cross-section.  This  is  consistent  with  the  results  re¬ 
ported  by  Forsberg  (Reference  6) , 

It  was  shown  in  Section  III  for  the  axi  symmetric  case  (n  =  0)  that  the  torsional  displacement 
uncouples  from  the  radial  and  axial  displacements..  One  of  the  eigenvalues  corresponds  to 
pure  torsional  vibration  while  the  other  two  involve  coupled  radial  and  axial  motion. 


A  simple  expression  for  the  torsional  vibration  can  be  obtained  using  Equations  115  and 
103.. 


A  (  X  ,  n 


X  (!-*)(  l  +  3k) 


After  solving  uiq  and  substituting  y.  Equation  119  reduces  to 


For  the  freely-supported  cylindrical  shell  a  boundary  condition  of  zero  tangential  displacement 
was  assumed  at  both  ends  of  the  shell  which  corresponds  to  clamped  ends  for  torsion  of  a 
hollow  shaft  or  tube. 


The  lowest  eigenvalue  for  the  axisymmetnc  modes  is  associated  with  radial  deformation 

(C>A)  for  a  value  of-^-<  2.  The  lowest  eigenvalue  is  associated  with  torsional  vibration  for 

ma  /  / 

values  of  the  parameter  >2.  For  values  of  —  greater  than  a  value  approximately  equal 
to  7 r,  the  largest  eigenvalue  corresponds  to  radial  deformation  (C>A).;  The  lowest  frequency 
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of  axisymmetric  vibration  for  a  freely-supported  cylinder  witli  -JL- >n  corresponds  to  pure 

m<A 

ton.^onal  motion.  1'he  highest  frequency  corresponds  to  radial  motion  (Figure  19).  This  result 
was  reported  by  Forsberg  in  Reference  6. 

For  small  values  of  -4*  and  low  values  of  n  the  second  eigenvalue  yields  amplitude  ratios 

(B>A,  C)  while  modes  having  larger  values  of  —  and  n  have  amplitude  ratios  (A>B,  C).  Thus 

m£i 

the  second  eigenvalue  corresponds  to  torsional  vibration  for  some  modes  and  axial  vibration 
for  other  modes,  as  indicated  in  Figure  17.  In  Figure  18  a  similar  result  is  obtained  for  the 
third  eigenvalue. 

It  can  be  seen  from  the  figures  that  for  the  freely-supported  cylinder  under  consideration, 
the  lowest  vibration  frequency  is  associated  with  radial  vibration  for  non-axisymmetric  modes 
except  for  long  cylinders  with  one  circumferential  wave..  Among  the  various  radial  vibration 
frequencies  the  lowest  value  for  the  unloaded  shell  is  obtained  for  one  axial  half-wave  (m  =  1). 

This  leads  to  the  largest  value  of  — .  The  number  of  circumferential  waves  associated  with 

b  ma 

the  lowest  radial  vibration  frequency  is  a  function  of  the  shell  geometry  or  length  to  radius 

0  cl 

ratio  —  and  radius  to  thickness  ratio—.  The  larger  the  ratio  of  radius  to  thickness  the  greater 
£l  n 

the  number  of  circumferential  waves  in  the  minimum  frequency  Vibration  mode. 

For  non-axisymmetric  vibration  frequencies  corresponding  to  the  second  and  third 
eigenvalues,  the  lowest  frequencies  are  obtained  for  modes  with  one  axial  half-wave  as  in  the 
case  of  radial  vibration,  and  increase  as  the  number  of  half-waves  increase.  The  number  of 
circumferential  waves  associated  with  the  lowest  frequency,  however,  for  a  given  number  of 
axial  half-wa~-es  is  one  and  the  frequency  increases  as  the  number  of  circumferential  waves 
increases.  Thus  for  non-axisymmetric  vibration  corresponding  to  the  second  and  third 
eigenvalues,  the  lowest  frequencies  are  those  for  one  circumferential  wave,  and  one  axial 
half-wave,  aud  may  involve  either  motion  which  is  axial  or  torsional  depending  on  the  vt  lue 

of  fer  a  given  shell.  For  cylinders  with  a  value  of^<  2  the  second  eigenvalue  corre¬ 
sponds  to  torsional  vibration  and  for  cylinders  with— >2  to  axial  vibration  (Figure  17). 

If  the  frequency  is  assumed  to  be  zero,  the  eigenvalues  for  any  mode  are  the  values  of  the 
buckling  parameter  for  that  mode.  From  these  eigenvalues  the  axial  loading  for  buckling 

given  the  symbol  q^Q,  can  be  determined  for  the  specified  mode 

A  (  X  ,n  )  =  q  \Z  (123) 

The  three  eigenvalues  A  (X,  n)  associated  with  each  deflection  mode  correspond  to  buckling 

deformations  which  are  radial,  axial,  or  tangential.  The  lowest  eigenvalue  for  non- 
axisymmetric  modes  is  associated  with  buckling  which  involves  radial  deformation  except 
for  the  modes  with  one  circumferential  wave. 


Since  the  eigenvalues  for  a  given  mode  shape  can  be  obtained  without  specifying  a  value  for 
q^  the  values  A (X,  n)  in  Equation  123  for  the  buckling  loads  are  the  same  as  those  for  the 

circular  frequency  of  an  unloaded  cylinder,  Equation  118, or  the  vibration  of  a  cylinder  loaded 
below  the  buckling  value.  Equation  117.  In  addition  the  amplitude  ratios  A/C  and  B/C  are  the 
same  for  the  various  displacements. 
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Both  sides  of  Equation  117  can  be  divided  by  A(X,  n)  to  obtain 


or 


2  2  ,  \Z 

y  u)  4-  q  2  \ 

A(X ,  n  ) 

2  2  .2 
r  W  q2  X 

A  (  X,  n)  ‘  A ( X, n) 


(  124) 


(125) 


After  substituting  forA(X, 
result  is  obtained: 


n)  in  Equation  125  from  Equations  118  or  123  the  following 


(126) 


or 


(  127) 


Equation  127  shows  for  a  given  cylinder,  and  for  the  values  of  A(X,  n)  associated  with  a 
given  deflection  mode,  that  the  square  of  the  vibration  frequency  varies  linearly  ,vith  the  axial 
load,  or  axial  load  parameter  q~.  For  positive  values  of 


q g  (compressive  loading*  the  square 


of  the  vibration  frequency  uj  varies  linearly  from  its  value  for  no  axial  loading,  u,o,  to  zero 
at  the  buckling  load  =  q^Q.  A  compressive  loading  beyond  this  point  would  lead  to  imaginary 


frequencies  and  indicate  that  the  cylinder  was  statically  unstable  for  that  deformation  mode 
and  loading.  Negative  values  of  correspond  to  a  tensile  loading  and  increase  the  natural 


vibration  frequencies  for  a  given  mode. 


The  relationship  between  the  vibration  frequencies  and  axial  load  given  by  Equation  127 
was  obtainable  because  the  deformation  mode  shape  for  natural  vibrations  and  buckling  are 
the  same  for  the  particular  boundary  conditions  considered.  The  frequency  and  axial  load 
parameters  were  combined  into  a  single  eigenvalue  expression. 


The  lowest  buckling  mode,  or  minimum  value  of  qgQ  is  generally  not  associated  vith  a  mode 

shape  having  one  axial  half-wave  as  is  the  case  for  the  minimum  vibration  frequency  for  no 
axial  load.  It  is  shown  in  Equation  123  that  the  minimum  buckling  load  q^Q  is  associated  with 

the  minimum  value  of  —  - .  For  a  given  number  of  axial  half-waves,  m,  the  lowest  vibration 


frequency  and  the  lowest  buckling  load  have  the  same  number  of  circumferential  waves,  n. 


The  buckling  values  q^  for  axial  compression  are  obtained  from  Equation  123  after  the 

eigenvalues  have  been  obtained  from  Equation  108  for  non-axisymmatric  modes  (n  ^  o)  and 
Equations  114  and  115  for  axisymmetric  modes.  The  buckling  results  obtained  are  the  same 
as  those  of  Flugge  (Reference  4)  since  the  equations  are  identical.  Sime  only  the  buckling 
envelope  is  given  in  Reference  4  for  various  radius  to  thickness  ratios,  the  complete  buckling 
curves  are  plotted  in  Figures  20  through  36  for  use  with  Equations  126  and  127. 

The  various  eigenvalues  were  obtained  from  solution  of  the  cubic  equation  except  where 
the  linear  solution  gave  greater  accuracy  for  the  lowest  root. 
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Tne  buckling  results  for  axial  compression  are  plotted  in  Figure  20  for  the  axisymmetric 
modes.  The  results  for  buckling  involving  radial  deformation  are  given  in  Figures  21  through 
34  for  radius  to  thickness  ratios  from  20  to  5000..  Figure  35  is  a  plot  of  the  buckling  values 
corresponding  to  the  second  eigenvalue  and  Figure  36  is  a  plct  of  the  buckling  values  corre¬ 
sponding  to  the  third  eigenvalue.  As  in  the  case  of  free  vibration  each  plot  for  the  second  and 
third  eigenvalue  contains  both  torsional  deformation  and  axial  deformation.  A  single  figure 
for  each  case  covers  the  range  of  radius  to  thickness  ratios  from  20  to  5000  for  these  higher 
buckling  values,  indicating  ;hai  bending  jffects  are  unimportant  for  these  modes. 


Since  buckling  modes  with  one  circumferential  wave  have  the  lowest  eigenvalue,  as  in  the 
vibration  case,  the  amplitude  ratio  for  radial  and  tangential  displacements  approaches  -1. 

These  displacements  are  larger  than  the  axial  displacement  for  values  of >3.5.  These 

buckling  modes  are  similar  to  column  buckling  with  little  deviation  from  the  circular  cross 
section.  In  Reference  4  Flugge  shows  that  the  buckling  values  for  this  mode  approach  the 

p 

Euler  column  values  asymptotically  for  large  values  of  — . 

For  the  case  of  axisymmetric  buckling  (n  =  0)  due  to  an  axial  load.  Figure  20  shows  that 
0 

for  values  of  2  the  lowest  buckling  value  is  associated  with  radial  deformation.  If  the 

values  of^j>2  the  lowest  axisymmetric  buckling  value  involves  torsional  deformation.  As 

in  the  vibration  case,  the  torsional  displacement  uncouples  from  the  radial  and  axial  displace¬ 
ment  for  the  axisymmetric  modes.  The  expression  for  this  pure  torsional  buckling  due  to 
axial  load  can  be  written. 


q 


20 


£ 

D 


(i-  u) 
2 


(  I  +  3k  ) 


For  values  of  k«l  this  can  be  approximated 


or 


q  -  s 

m20  d  2 


2(l+i/) 


(120) 


(129) 

(130) 


in  terms  of  the  axial  load  per  unit  length  on  the  cylinder  boundary.  In  Equation  129  it  is  shown 
that  the  value  of  q^Q  for  torsional  buckling  due  to  an  axial  load  is  a  constant  for  a  given  cyl¬ 
inder.  This  also  is  shown  in  Figure  20. 


While  values  of  the  buckling  parameter  q^Q  for  various  modes  corresponding  to  the  second 

and  third  eigenvalue  are  unimportant  from  a  stability  standpoint  they  are  plotted  in  figures  35 
and  36  for  use  with  Equation  126  or  127.  As  in  the  vibration  case  the  minimum  buckling  values 
for  the  second  and  third  eigenvalues,  for  the  non-axisymmetric  modes,  occur  for  one  circum¬ 
ferential  wave  and  sn  ail  values  of  —  or  a  large  number  of  axial  half-waves. 

The  vibration  frequencies  for  freely-supported  cylinders  with  an  axial  load  can  be  deter¬ 
mined  by  using  Figures  20  through  36  and  Equation  126  or  127..  These  equations  are  valid 
for  any  of  the  three  eigenvalues  for  a  given  detormation  mode.  For  a  given  cylinder  it  would 
be  necessary  to  select  a  mode  of  interest  by  assigning  integer  values  to  n  and  m,  and  calcu¬ 
lating  a  value  for  After  determining  whether  radial,  torsional,  or  axial  vibration  is  of 
interest  the  values  of  and  can  be  determined  Irom  the  appropriate  charts  for  that 
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particular  n  and  The  value  of  the  parameter  for  the  applied  loading  can  be  calculated 

and  the  new  frequency  u)  or  the  frequency  ratio  can  be  calculated  from  Equation  126  or 

% 

127,  If  the  shell  is  subjected  to  an  axial  tension  loading,  a  negative  value  for  q  is  used  in 
the  formulas  which  yields  an  increase  in  frequency. 

It  can  be  seen  that  the  greatest  change  in  frequency  occurs  in  those  modes  having  the  lowest 
buckling  values  q^ 
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AXIAL  MODE  SHAPES 
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Figure  2.  Deflection  Modes 
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Figure  3.  Radial  Vibration  Frequencies.  -jp  20.  No  external  Pressure  (c^-  0) 
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Radial  Vibration  Frequencies.  20.  No  external  Pressure  (q^  =  0} 
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Figure  4. 
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Figure  8. 


Radial  Vibration  Frequencies. 
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100.  No  External  Pressure  (q^  :  0) 
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Figure  J,  Radial  Vibration  Frequencies,  -jp  500.  No  External  Pressure  (q  ; 
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NUMBER  OF  CIRCUMFERENTIAL  WAVES  n 


Figure  10.  Radial  Vibration  Frequencies.  -j--500.  No  External  Pressure  (q^  -  0) 
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Figure  12.  Radial  Vibration  Frt*quencies.  -j^lOOO.  No  External  Pressure  (qj  :  0) 
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Figure  13.  Radial  Vibration  Frequencies.  2000.  No  External  Pressure  (q 
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Figure  14.  Radial  Vibration  Frequencies.  ■^-I2000.  No  External  Pressure  (q,  =  0) 
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Figure  15. 
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Figure  16.  Radial  Vibration  Frequencies. -^-r5000.  No  External  Pressure  (q^  0) 
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Second  Vibration  Frequencies,  -jp- 20-5000.  No  External  Pressure  (q^-  0) 
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EIGENVALUE 


AFFDL-TR- 67-28 


Figure  18.  Third  Vibration  Frequencies.  — -  2<>-5000.;  Wo  External  Pressure  (q^  =  0) 
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Figure  Id.  Axisymmetric  Vibration  Frequencies.  -j*-1  20-5000.  No  External  Pressure 
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Figure  20. 
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Figure  21.  Buckling  Values  for  Axial  Compression.  ~?~r20.  No  External  Pressure 
(qj  :  o) 
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Figure  22.  Buc ing  Values  for  Axial  Compression.  20.  No  External  Pressure 
(qj  s  V 
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j  23  Buckling  Values  for  Axial  Compression.— ■=  50.  No  External  Pressure 


BUCKLING  LOAD  PARAMETER 


BUCKLING  LOAD  PARAMETER 


BUCKLING  LOAD  PARAMETER 


AFFDL-TR-67-2c< 


NUMBER  OF  CIRCUMFERENTIAL  WAVES  n 


Figure  26.  Buckling  Values  for  Axial  Compression.  ^-=100.  No  External  Pressure 
(Qj  =  0) 
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Figure  27.  Buckling  Values  for  Axial  Compression.  -r-~  500.  No  External  Pressure 

(q,  =  o)  b 
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Figure  28.  Buckling  Values  for  Axial  Compression.  -^-=500.  No  External  Pressure 
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Figure  31.  Buckling  Values  for  Axial  Compression.  —=2000.  No  External  Pressure 
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Figure  32.  Buckling  Values  for  Axial  Compression.  -r~;2000.  No  External  Pressure 
(qj  --  0) 
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Figure  34. 


Buckling  Values  Tor  Axial  Compression.  -r-r  5000.  No  External  Pressure 
(qj  --  0) 
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Figure  35.  Second  Buckling  Value  for  Axial  Compression.  -g-  =  20-5000.  No  External 
Pressure  (q,  =  0) 
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Figure  3S.  Third  Buckling  Value  for  Axial  Compression. -^—=20-5000.  No  External 
Pressure  (q  -  0)  n 
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SECTION  V 

BUCKLING  RESULTS  FOR  EXTERNAL  PRESSURE 


The  buckling  values  f~r  a  freely-supported  cylinder  subjected  to  only  external  pressure  can 
be  obtained  by  setting  the  parameters  q ^  and  oi  equal  to  zero  in  the  various  terms  of  the  char¬ 
acteristic  determinant  and  separating  the  terms  containing  q^  from  the  expressions  D  in 

Equation  98.  The  external  pressure  parameter  q.  is  given  the  symbol  q  ^  to  indicate  that  it 
is  the  buckling  value. 


‘  <Hu  -  9,0  »2  > 

H  12 

^ R.3  —  Qio  X ) 

\A  1 

H|2 

(H22“  9.0  "2  ] 

(H23~ql0  n) 

8 ! 

>  =0 

(131) 

-  ^13  ~  9io  X  ) 

<H«-q,Ofl> 

(H33-q,o  n  )_ 

(  c 

) 

Expressions  for  the  various  terms  H.j  can  be  obtained  by  subtracting  the  q^  terms  from 
the  expressions  for  in  Equations  100  through  105. 

For  a  non-trivial  solution,  the  determinant  of  the  3x3  matrix  above  must  be  zero.  After 
expanding  this  determinant  and  collecting  like  powers  of  q^  the  following  third  order  char¬ 
acteristic  equation  for  buckling  due  to  external  pressure  is  obtained:. 

[~"6  +  X*  n2  ♦  q*o  [nMH^H^+H^  )+2H|2nA  -2n2\H|3-\2H22 
-2n3H23-nZ  H„  ]  +  q,0  [-  n2  ( H„  «22  +  H„  H„  +  H22  H„)  +  n2  (  H  f3  +  +  H  J, ) 


-2H,2(n  H|3+\H23)+2XH13H22  +  2n  HI|H23]+  H„  HZ2H33+2HI2  Hl3  H„ 
H2  2  H  33  ~H23  H„  --0 


(132) 


Equation  132  is  valid  only  tor  non-axisymmetric  modes  (n>o).  As  in  the  case  of  vibration  and 
buckling  due  to  an  axial  load,  the  characteristic  equation  for  buckling  due  to  external  pressure 
is  a  cubic  equation,  and  three  buckling  values,  q^Q,  can  be  obtained.  The  expressions  for  the 

amplitude  ratios  for  the  various  buckling  values  are 


_A 

C 


(H22  ql0n  )  (H|3  q ,0 X  ) 


~  <H23  “<»ion  > 


M2 
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As  for  both  vibration  and  buckling  due  to  axial  'oad,  the  non-axi symmetric  modes  involve  a 
a  coupling  between  the  three  displacements  u,  v,  and  w. 

For  some  modes,  the  solution  of  Equation  132  yielcs  one  negative  and  two  positive  buckling 
values.  The  negative  value  for  indicates  bucklingdue  to  internal  pressure,  since  a  positive 

value  for  q1Q  corresponds  to  external  pressure.  In  those  cases  where  a  negative  buckling 

value  was  obtained  for  q^  there  were  two  buckling  values  for  the  given  mode  which  had 

primarily  radial  deformation.  One  of  these  was  associated  with  the  negative  value  of  and 

the  other  with  the  lowest  positive  value  for  q^.  The  only  exception  is  for  modes  having  one 

circumferential  wave,  where  for  long  shells  the  radial  and  tangential  displacements  were 
approximately  equal  and  both  were  greater  than  the  axial  deformation  (C*B>A),  which  would 
represent  “beam  modes”  for  the  cylinder.  In  the  cases  where  three  positive  values  were 
obtained  for  q10  one  corresponded  to  radial  deformation  (C>A,  B)  one  to  tangential  deforma¬ 
tion  (B>A,  C)  and  one  to  axial  deformation  (A>B,  C). 


The  lowest  positive  buckling  values  for  q^  (external  pressure)  for  various  modes  are 

plotted  in  Figures  37  to  50  for  various  values  of  radius  to  thickness  from  20  to  5000.  There 
are  two  figures  for  each  radius  to  thickness  ratio.  In  one  figure  qn  ft  is  plotted  versus  the 

Q  iU 

axial  wavelength  parameter  for  various  values  of  n,  and  in  the  other  qin  is  plotted  versus 

the  number  of  circumferential  waves,  n,  for  various  values  of  - —  .A  value  of  the  lowest  buck- 

ma 

ling  load  for  the  various  modes  was  obtained  from  the  cubic  characteristic  equation.  Equa¬ 
tion  132,  and  also  the  linear  form  of  this  equation,  and  the  more  accurate  result  plotted  in  the 
figures,  "Hus  was  determined  by  substituting  each  result  back  into  Equation  132  and  checking 
the  residual.  For  small  values  of  q^Q  associated  with  radius  to  thickness  ratios  from  2000  to 

5000  and  values  of  —  >5,  the  linear  equation  generally  gave  greater  accuracy  for  the  lowest 

root.  This  is  because  of  the  numerical  accuracy  which  was  obtained  from  the  digital  computer 
program  used  to  obtain  the  roots  of  the  cubic  equation.  As  for  both  vibration  and  buckling  due 
to  axial  load,  the  lowest  values  of  the  character  stic  equation  were  generally  associated  with 
deformation  which  was  primarily  radial  and  the  results  were  dependent  on  the  radius  to 

thickness  ratio  ^  indicating  that  bending  effects  are  important  for  these  cases. 


The  modes  having  one  circumferential  wave  do  not  appear  in  Figures  37  through  50  because 
the  buckling  values  for  external  pressure  for  these  modes  are  higher  than  for  modes  with  two 
or  more  circumferential  waves.  The  lowest  buckling  value  for  external  pressure  decreases 

as  the  parameter increases,  indicating  that  for  a  given  cylinder  the  lowest  buckling  value, 

q^g,  is  associated  with  a  mode  having  one  axial  half-wave,  since  this  leads  to  the  largest 

value  of — ^  .  For  the  large  values  of  the  lowest  buckling  values  for  external  pressure 

are  associa  ed  with  two  circumferential  waves  (n  =  2).  As  the  value  of  the  parameter-—-  de- 

ma 

creases  the  cumber  of  circumferential  waves  associated  with  the  lowest  buckling  value  in¬ 
creases.  These  results  were  obtained  by  Flugge  (Reference  4),  using  only  the  linear  terms  of 
the  characteristic  equation.  Flugge  points  out  Uiat  the  linear  form  of  the  characteristic  equa- 

h2 

tion  n  =  2  and  X=  0,  where  X  =  gives  a  value  of  3k  for  qin  where  k  = - 5-  indicating 

*  1U  12a 


that  the  curve  for  n  =  2  is  asymptotic  to  3k  for  large  values  of 


ma 


57 


AFFDL-TR-67-28 


The  buckling  results  for  external  pressure  for  modes  having  one  circumferential  wave  are 

given  in  Figure  51  for  values  of  from  20  to  5000.  Two  buckling  values  were  obtained  for 

external  pressure  and  one  negative  buckling  value  for  q1Q  which  represents  internal  pressure. 

For  values  of  —  <  tt  there  are  two  buckling  modes  which  involve  radial  deformation,  with  the 
ma 

value  for  external  pressuxe  lower  than  the  value  for  internal  pressure.  There  is  also  a  higher 
buckling  value  for  external  pressure  which  involves  torsional  deformation.  The  buckling  values 

appear  to  be  a  function  of  the  ratio  of  radius  to  thickness  for-— ^<1,  and  are  essentially  the 

same  for  values  of  from  20  to  5000  for-^-  >  1. 

h  ma 

Flugge  (Reference  4),  in  working  with  the  linear  form  of  the  characteristic  equation,  points 
out  that  the  buckling  load  qin  for  this  case  changes  from  positive  to  negative  at  —  =  if,  with 

A  v*  HAS 

the  negative  values  associated  with  — -  >7 r.  In  Figure  51,  where  all  three  buckling  values  are 

nia. 

shown,  it  can  be  seen  that  the  buckling  curve  for  internal  pressure  crosses  the  buckling  curve 

for  external  pressure  at  -JL  =  tt.  Thus  the  lowest  buckling  load,  q,n,  occurs  for  external 
0  ma  a  6  10  Q 

pressure  for  —  <  7r  and  internal  pressure  for-=->  7 r.  For  values  of  —  >  ir  the  buckling 
ma  ma  ma  0 

results  for  internal  pressure  have  an  amplitude  ratio  B/C  which  is  approximately  -1,  indicating 
that  the  amplitude  of  the  radial  and  tangential  displacements  are  approximately  equal  and  both 
are  larger  than  the  axial  deformation.  This  represents  deflection  with  little  distortion  of  the 
cylinder  cross  section  and  would  represent  “beam  modes’*  for  the  cylinder. 

For  values  of -^>4  the  lowest  buckling  value  for  external  pressure  represents  axial  de¬ 
formation.  Thus,  for  values  of  ^~>4  the  buckling  value  for  internal  pressure  was  a  “beam 

mode”,  while  the  lowest  buckling  value  for  external  pressure  involved  axial  deformation,  and 
the  second  buckling  value  for  external  pressure,  torsional  deformation. 


For  the  case  of  axisymmetric  buckling  (n  =  0  modes)  due  to  external  pressure,  the  terms 
H^2  &nd  in  Equation  131  are  zero,  and  the  tangential  displacement  v  can  be  uncoupled 

from  the  axial  and  radial  displacements  as  for  both  vibration  and  buckling  due  to  axial  load. 
Equation  98  can  be  written  as  two  equations  for  these  modes. 

^13  ”  ^  9,0  H33 

H22  B  =0  (136) 


The  uncoupled  equation  for  torsional  deformation.  Equation  136  has  only  the  trivial  solution 
B  =  0,  indicating  that  axisymmetric  torsional  buckling  does  not  exist  for  a  freely-supported 
cylinder  subjected  to  external  pressure. 


The  determinant  of  the  2x2  matrix  in  Equation  135  is  expanded  to  form  the  quadratic 
characteristic  equation  for  the  axisymmetric  case. 


<0  (Xl)  -<,0(2XH,3)+  <3 


-HUH33  10 


( 1 37) 
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The  amplitude  ratio  involves  only  coupled  radial  and  axial  displacement  for  this  case. 


A  f  H  (  3 

_  — 


(138) 


The  huckling  values  for  external  pressure  obtained  from  Equation  137  are  plotted  in  Figure  52. 
One  positive  value  for  q1Q  and  one  negative  value  for  q1Q  were  obtained  for  range  of  radius  to 

thickness  ratios  from  20  to  5000  and  the  range  of  the  axial  wavelength  parameter,  from 

.2  to  15.  The  lower  buckling  value  is  for  external  pressure.  For  small  values  of  k,  (k<<l), 

.2 

where  k  = - g-,  and  for  X  =  0,  where  X  =  “j>i  Equation  137  yields  the  following  values  for 


12a 


no' 


±i 


039) 


This  gives  a  buckling  value  for  external  pressure  of  q1Q  =  1-  v,  and  a  buckling  value  of 
a  =  -1-v,  which  represents  internal  pressure,  for  I^ng  thin  shells.  If  Poisson’s  ratio.y,  were 
zero,  the  buckling  value  for  internal  and  external  pressure  would  be  1.  In  Figure  52  it  is 

shown  that  these  simple  buckling  values  appear  to  be  valid  for  values  off- from  500  to  5000 

/  /  an 
for  the  full  range  of  —  examined  (,2<  — -  <  15),  and  for  values  of-r  from  20  to  100  for 
**  ma  '  ma  '  h 

—  >  1. 
ma 

The  buckling  values  for  internal  pressure  for  all  modes  are  shown  in  Figure  53.  These 

£ 

buckling  values  are  associated  with  radial  deformation,  except  for  n  ~  1  with  —  >  v,  which 
represents  the  “beam  mode”  as  previously  discussed.  The  modes  with  two  circumferential 

p 

waves  yield  a  negative  value  for  q^  for  1,  which  appears  to  get  very  large  as  ap¬ 


proaches  1.  For  a  value  of^~  >1,  Equation  132  gave  only  positive  buckling  values  for  q^. 
Similarly  for  three  circumferential  waves,  a  negative  buckling  value  could  be  obtained  for 

Q  g  Q 

—  <  .4  and  all  values  of  —  examined,  but  at  >,4  only  positive  values  are  obtained  for  q^. 


ma 


ma 


Although  it  is  not  shown  on  the  figure,  a  large  negative  value  for  q.n  was  obtained  for  four 

£  /  lu 

circumferential  waves  for  —  =  .2,  while  for  all  values  of  —  ^  .25  only  positive  buckling 

ma  ma 


values 


*10 


were  obtained. 


The  second  positive  buckling  value  for  external  pressure  is  plotted  in  Figure  54,  and  the 
third  positive  buckling  value  is  plotted  in  Figure  55.  As  for  both  vibration  and  buckling  due  to 
axial  load,  the  second  buckling  value  for  external  pressure  involves  torsional  deformation  for 
some  modes  and  axial  deformation  for  the  rest.  A  similar  result  is  obtained  for  the  third 
buckling  value  for  external  pressure.  The  second  and  third  positive  buckling  values  were 
essentially  the  same  for  all  radius  thickness  ratios  from  20  to  5000,  indicating  the  predomi¬ 
nance  of  membrane  effects.  In  Figure  55,  which  gives  the  third  buckling  value,  the  curves  for 

n  =  2,  3,  and  4  were  stopped  at-^-  =1,  .4,  and  .25  respectively,  since  a  third  positive  value 


ma 


for  q,n  did  not  exist  for  the  next  increment-^-  below  these  values, 
10  ma 
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Figure  37.  Buckling  Values  for  External  Pressure.-^  20.  No  Axial  Load  (q.f  =  0) 
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Figure  38.  Buckling  Values  for  External  Pressure. -g-- 20.  No  Axial  Load  - 
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Figure  39.  Buckling  Values  for  External  Pressure. -^-=50.  No  Axial  Load  (q^  =  0) 
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Figure  40.  Buckling  Values  for  External  Pressure.  50.  No  Axial  Load  (q„  =  0) 
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Figure  41.  Buckling  Values  for  External  Pressure,  y- 100.  No  Axial  Load  (q£  =  0) 
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Figure  43.  Buckling  Values  for  External  Pressure,  -g-1 500.  No  Axial  Load  (q9=  0) 
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Figure  44.  Buckling  Values  for  External  Pressure.  —-'500.  No  Axial  Load  (q^  r  0) 
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Figure  45.  Buckling  Values  for  External  Pressure.  ^  :  iOOO,  No  Axial  Load  =  0) 
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Figure  46.  Buckling  Values  for  External  Pressure.  •— =  1000.  No  Axial  Load  (q2  :  0) 
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Figure  47.  Buckling  Values  for  External  Pressure.  -|-=2000.  No  Axial  Load  (q2  =  0) 


70 


BUCKLING  LOAD  PARAMETER  q 


AFFDL-TR-67-28 


71 


BUCKLING  LOAD  PARAMETER  q 


AFFDL-TR-67-28 


Figure  49.  Buckling  Values  for  External  Pressure.  — -5000. No  Axial  Load  (a,  =0) 
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Figure  50.  Buckling  Values  for  External  Pressure.  5000.  No  Axial  Load  (q2  =  0) 
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Figure  51.  Buckling  Values  for  One  Circumferential  Wave,  Internal  and  External 
Pressure.  20-5000.  No  Axial  Load  (q9  =  0) 
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Figure  52.  Axisymmetric  Buckling,  Internal,  and  External  Pressure,  -r--  20-5000.  No 
Axial  Load  (q„  1  0) 
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Figure  54. 
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<q2  :  0} 


77 


ING  LOAD  PARAMETER  g,0 


AFFDL-TR-67-28 


(q2  -  0) 


78 


AFFDL-TR-67-2S 


SECTION  VI 

THE  VIBRATION  OF  FREELY-SUPPORTED  CYLINDERS 
SUBJECTED  TO  AXIAL  LOAD  AND  EXTERNAL  PRESSURE 


It  was  shown  In  Section  IV  that  the  axial  load  parameter,  q 2  can  be  combined  with  the  fre¬ 
quency  parameter,  oi,  to  form  a  single  eigenvalue  expression,  while  the  internal  pressure 
parameter  appears  in  the  coefficients  of  the  resulting  characteristic  equation.  An  exami¬ 
nation  of  the  expressions  P.^,  given  in  Equations  100  through  105,  which  appear  in  Equation  98, 

2 

reveals  that  the  term  -q^n  appears  on  the  main  diagonal  of  that  equation.  These  terms  can 


be  moved  to  the  right-hand  side  of  Equation  98  and  combined  with  the  eigenvalue  expression 
as  was  done  with  the  parameter  q^.  The  factor  q^  also  appears  in  the  terms  and  D23, 


indicating  that  the  coefficients  of  the  resulting  characteristic  equation  would  still  be  a  function 


m 

'13 


3  external 

pressure  parameter  q^.  With 

these 

»: 

um 

HI2 

(H.'J-q,  X) 

t^m| 

12 

"22 

(Hjj  -  q,  n) 

jBm 

,  -  q,  X) 

(H»  -q(  ?> ) 

H  33 

(cm 

‘m 


kB„?s(yV+<,*+V8>  ?|Bfn 

C 


(140) 


m 


where  the  terms  H..  are  given  by  subtracting  the  q1  terms  from  the  expressions  for  D..  (Equa¬ 
tions  100  to  105).  11  L  3 


If  the  value  of  q^^  is  small  in  the  expressions  for  H13  and  H,,^,  the  parameter  q^  can  be 

dropped  from  the  left-hand  side  of  Equation  140  without  affecting  results  significantly.  With 
this  simplification,  a  new  eigenvalue  expression  A(X,  n)  would  result 

A  (  X  ,  n  )  =  (  yZ  of  4-  q2  X2  +  q(  n2  )  (141) 

The  eigenvalues  A (X,  n)  would  be  the  same  as  A(X,  n)  of  Section  IV  because  the  coefficients 
in  Equation  140  are  the  same  as  the  with  qj  =  0. 

A  buckling  value,  q^,  can  be  obtained  from  Equation  141  bv  setting  ut  and  q^  equal  to  zero. 

A  (X  ,n)  =  A(  X  ,n  )  =q(Qn2  (142) 


q 


10 


A  (  X.n  ) 

2 

n 


(143) 


If  the  approximation  discussed  in  Section  IV  is  assumed  to  be  valid,  an  expression  similar  to 
Equations  124  and  125  which  include  the  effects  of  both  axial  load  and  internal  pressurization 
can  be  obtained. 


( 144) 
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or 


2 


UJ 


!  145) 


Equations  144  and  145  indicate  that  the  square  of  the  vibration  frequency  for  a  freely- 
supported  cylinder  varies  linearly  with  both  axial  load  and  internal  pressure  for  a  given 
deformation  mode.  For  buckling  due  to  combined  loads.  Equation  145  can  be  written 


I  = 


( 146) 


These  results  are  shown  in  Figure  56.  The  line  (— )=  0  represents  buckling  of  the  cylinder. 

wo 

Thus  the  coordinates  of  any  point  on  this  line  represent  combinations  of  axial  loading  and 
external  pressure  necessary  to  buckle  the  cylinder  in  a  given  mode.  Similarly  the  line 

~  0.5  would  represent  the  combinations  of  axial  load  and  external  pressure  which  reduce 

the  square  of  the  vibration  frequency  to  a  value  of  one  half  that  for  the  unloaded  cylinder.  The 

line(~ -)  =  1  represents  combinations  of  axial  load  and  external  pressure  for  which  the  vibra¬ 
to 

tion  frequency  remains  unchanged. 


The  results  obtained  for  buckling  due  to  external  pressure  in  Section  V  indicate  that  for 
some  modes  negative  buckling  values  are  obtained,  which  correspond  to  buckling  due  to  internal 
pressure.  With  Equation  143,  however,  only  positive  buckling  values  could  be  obtained  since 
all  the  eigenvalues  A(X,  n)  are  positive.  Thus,  buckling  values  obtained  using  Equation  143 
could  lead  to  significant  errors  for  some  modes  If  the  actual  results  for  q10  obtained  from 


the  buckling  equations  of  Section  V  are  used  in  Equations  144  or  145  rather  than  the  values 
given  by  Equation  143,  the  problems  arising  from  the  negative  buckling  values  would  be  al¬ 
leviated.  The  use  of  the  true  values  for  q^  obtained  from  the  buckling  equations,  essentially 

^li 

forces  Equations  144  and  145  to  be  correct  at—  =  1.  In  addition,  use  of  the  actual  buckling 

q10 

results  for  q^Q  was  found  to  exterd  the  range  ol  validity  of  Equations  144  and  145  considerably, 


that  is,  the  equations  were  found  to  be  valid  where  the  buckling  value  was  of  the  same  order 
of  magnitude  as  v. 


It  is  proposed  that  Equations  144,  145  cr  Figure  56  be  used  together  with  Figures  3  through 
55  to  obtain  values  for  the  vibration  frequencies  or  buckling  values  for  freely-supported  cyl¬ 
inders  subjected  to  axial  loading  and  external  pressure.  For  a  given  freely-supported  cylinder, 

D 

the  ratio  of  length  to  radius, and  radius  to  thickness,  t  would  be  fixed.  A  particular  mode  of 

interest  could  be  selected  by  assigning  integer  values  to  the  number  of  circumferential  waves, 

n,  and  the  number  of  axial  half-waves,  m,  in  the  parameter  It  would  then  be  necessary  to 

establish  which  of  the  eigenvalues  associated  with  the  particular  mode  shape  is  to  be  examined; 
that  is,  whether  radial,  torsional,  axial,  or  in  some  cases,  the  “beam”  modes  are  of  interest. 
The  values  of  q^,  and  q^Q  for  the  particular  mode  of  interest  can  be  obtained  from  the 

appropriate  figure.  If  a  given  axial  loading  or  external  pressure  is  specified,  a  value  for  the 
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parameter  q 2  or  qj  can  be  calculated.  With  this  information  the  vibration  frequency  or  reduc¬ 
tion  in  frequency  can  be  obtained  from  Equation  144, 145,  or  Figure  56.  For  combined  buck¬ 
ling  results  either  q^  or  q^  could  be  taken  as  the  unknown  to  be  determined.  The  greatest 

change  in  vibration  frequency  occurs  for  those  modes  having  the  lowest  buckling  values  *10 
and  which  are  usually  those  involving  primarily  radial  deformation. 

No  simplification  or  approximations,  other  than  those  inherent  in  the  original  differential 
equations  or  in  the  assumption  of  the  form  of  the  solution,  are  associated  with  the  axial  load 
effects,  discussed  in  Section  IV.  The  simplifications  made  in  arriving  at  Equations  144  to  146 
or  Figure  56  involve  only  the  external  pressure  effects.  The  range  of  validity,  or  errors  re¬ 
sulting  from  use  of  the  procedures  previously  outlined,  were  determined  by  comparing  these 
results  with  those  obtained  by  using  Equation  98  or  108  with  all  of  the  q^  terms  retained.  To 

avoid  errors  arising  from  reading  the  figures,  the  digital  results,  used  in  plotting  the  figures, 
were  used  in  this  error  analysis.  Calculations  were  made  for  the  full  Tange  of  the  various 
parameters  involved. 


Buckling  calculations  were  made  for  values  of  pressurization  equal  to  *.3qjp,  .3q^Q,  .6q^ 
and  1.4q^Q  for  various  modes.  The  value  of  used  in  the  calculations  was  that  which  was 
smallest  in  absolute  value  for  that  mode.  In  some  cases  this  was  a  negative  value,  associated 

^2 

with  internal  pressure.  The  values  of —  for  buckling  associated  with  the  above  pressurization 

values  were  determined  by  the  two  methods.  A  comparison  of  these  results  was  made  for  0, 
1,  2,  5,  10,  and  20  circumferential  waves,  values  of  equal  to  .2,  .5,  1,  5,  and  10,  and  for 

EQ8. 

radius  to  thickness  ratios  from  20  to  5000.  Since  the  linear  relationship  between  the  square 
of  the  vibration  frequency  and  axial  load  was  established  without  any  approximations,  the 

^2 

relative  error  in  —  for  combined  buckling  would  be  the  same  as  for  the  square  of  the  vi- 
**20 
2 

bration  frequency  for  ladings  less  than  the  buckling  values. 

In  general  a  comparison  of  the  results  indicates  that  for  the  bulk  of  the  cases  examined  the 

% 

difference  in  the  values  obtained  for  —  were  less  than  1%  of  q_n.  The  largest  discrepancy 

q20 

occurred  for  the  axisymmetric  modes  (n  =  0',  and  for  modes  with  one  circumferential  wave. 

The  results  for  the  axisymmetric  modes  are  shown  in  Figure  57.  The  plotted  points  represent 

ql  ql 

the  correct  results.  The  two  methods  produce  the  same  results  at  — ~  =0  and  — —  =  1.  The 

a  1  qi°  qi° 

greatest  discrepancy  occurred  for  ~  =  20  and  — -  =  .2,  whichwould  represent  a  relatively  thick 

u  in& 

short  cylinder.  The  results  for  —  =  5  wore,  essentially  the  same  for  all  values  of  the  radius 
to  thickness  ratio  examined  and  produced  the  smallest  discrepancy. 

The  error  is  generally  less  ac  n  increases.  Figure  58  shows  a  comparison  of  the  results 

n 

for  modes  with  one  circumferential  wave.  Again  the  largest  discrepancy  occurs  for  -j-  -  20 

and  =  .2.  For  values  ot-^-  >  5  the  error  was  less  than  2.3%  of  q»ft  for  ail  values  of  radius 
ma  ma  ^20 

to  thickness. 
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The  results  for  two  circumferential  waves  are  shown  in  Figure  59.  For  values  of  — •  equal 

iUd 

to  .5  and  1,  the  errors  ranged  from  .59  to  6.8%  of  for  all  values  of  the  radius  to  thickness 

ratio.  For  values  of — ->  5  the  errors  were  much  less  than  1%  for  all  values  of  the  radius  to 
thickness  ratio.  ma 


a  / 

For  five  circumferential  waves  the  error  for  =  20  and—  =  .2was  1.7%  of  at-—  ~  ”.3* 

q-j 

92%  of  q_n  at  —  =  .3,  1.05%  of  q_n  at  =  .6,  and  2.39%  of  q~n  at  =  1.4.  and  were  less 
20  q1Q  2U  qin  m  q1ft 


no 


10 


Q  it 

less  than  1%  for  all  other  values  of--  and  v-,  and  in  most  cases  substantially  less  than  1%. 

ma  h 

For  modes  having  10  and  20  circumferential  waves  the  results  were  found  to  agree  to  well 

9  ^ 

within  1%  for  the  entire  range  of  ~  and  p 


82 


RATIO  OF  AXIAL  LOAO  TO  BUCKLING  VALUE 


AFFDL-TR-67-28 


Figure  56.  Reduction  in  Vibration  Frequencies  for  Combined  Effects  of  Axial  Load  and 
Pressurization 
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Figure  57.  Comparison  of  Combined  Buckling  Results  for  Axisymmetric  Modes,  (n  =  0) 
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Figure  58 .  Comparison^  Combined  BuckUng  Results  for  Modes  with  One  Circumferential 
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Figure  59.  Comparison  of  Combined  Buckling  Results  for  Modes  with  Two  Circumferential 
Waves  (n  =  2) 
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SECTION  VII 

RESULTS  AND  CONCLUSIONS 


Equations  are  formulated  which  can  be  used  to  determine  the  vibration  and  buckling  char¬ 
acteristics  of  cylindrical  shells  with  various  homogeneous  boundary  conditions  under  axial 
loading  and  a  constant  external  pressure.  A  method  of  solution  is  outlined  which  makes  use  of 
a  series  of  beam  f  motions  for  the  radial,  tangential,  and  axial  displacements,  which  satisfy 
the  appropriate  boundary  conditions  involving  these  displacements. 


Using  a  series  of  functions  known  to  be  an  exact  solution  to  the  problem  for  the  freely- 
supported  cylinder  the  equations  reduced  to  those  obtained  by  Flugge  for  the  buckling  of  a 
cylinder  subjected  to  axial  loads  and  a  constant  external  pressure,  if  the  frequency  is  assumed 
to  be  zero.  If  the  axial  load  and  external  pressure  parameters  are  assumed  to  be  zero,  the 
equations  reduce  to  a  form  similar  to  those  obtained  by  Arnold  and  Warburton  for  the  natural 
vibrations  of  a  freely-supported  cylinder. 


An  investigation  was  made  of  the  vibration  and  buckling  characteristics  for  the  freely- 
supported  cylinder  under  axial  loading  and  external  pressure,  since  an  exact  solution  was 
shown  to  exist  for  this  case.  The  results  discussed  in  the  remainder  of  this  Section  are 
associated  with  this  solution. 

The  characteristic  equations  for  natural  vibrations,  buckling  due  to  axial  load,  and  buckling 
due  to  external  pressure  were  solved  for  a  range  of  cylinder  raoius  to  thickness  ratios  from 

20  to  5000,  for  values  of  the  parameter  from  .2  to  15  and  for  all  circumferential  waves 

ma 

from  0  to  20.  The  vibration  and  buckling  results  are  presented  in  a  series  of  figures. 

For  the  freely-supported  cylinder  the  problem  reduces  to  the  solution  of  a  cubic  character¬ 
istic  equation,  except  for  the  axisymmetric  modes,  indicating  that  three  eigenvalues  exist  for 
each  mode.  A  particular  mode  involves  an  integer  number  of  circumferential  waves  and  axial 
half-waves.  These  eigenvalues  are  associated  with  various  amplitude  ratios  for  the  radial, 
axial,  and  tangential  displacements,  indicating  thai  these  displacements  are  coupled.  For 
vibration  or  buckling  under  axial  loading  three  positive  eigenvalues  are  obtained  for  a  given 
mode  and  are  associated  with  deformations  which  are  primarily  radial,  axial,  and  torsional. 
For  long  cylinders  having  one  circumferential  w.'.ve  the  amplitude  of  the  radial  and  tangential 
displacement  approaches  a  value  of  -1  for  one  of  the  eigenvalues.  This  represents  the  “beam” 
modes  for  the  cylinder,  where  the  cylinder  deflects  with  little  distortion  of  the  circular  cross- 
section. 

For  buckling  due  to  external  pressure  a  negative  eigenvalue  is  obtained  for  some  modes, 
which  corresponds  to  buckling  due  to  internal  pressure.  In  this  case  the  negative  eigenvalue 
and  one  of  the  positive  eigenvalues  for  that  mode  are  both  associated  with  deformation  which 
is  primarily  radial. 

For  the  axisymmetric  mode  (n  =  0),  the  tangential  displacement  uncouples  from  the  axial  and 
and  radial  displacements.  For  the  natural  vibration  case  and  buckling  due  to  axial  load,  the 
uncoupled  torsional  displacement  gives  rise  to  a  pure  torsional  vibration  or  buckling  mode. 
Two  eigenvalues  are  associated  with  the  coupled  radial  and  axial  displacements.  For  buckling 
due  to  only  external  pressure,  the  uncoupled  torsional  deformation  has  only  the  trivial  solu¬ 
tion,  and  only  the  two  eigenvalues  associated  with  coupled  radial  and  axial  displacement 
can  be  obtained. 
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The  eigenvalues  associated  with  radial  deformation  were  generally  dependent  on  the  radius 

to  thickness  ratio  -p  whereas  the  eigenvalues  associated  with  axial  or  torsional  deformation 

were  independent  of  the  radius  to  thickness  ratio  for  values  of  ---from  20  to  5000. 

h 

The  third  order  characteristic  equation,  the  resulting  eigenvalues,  and  amplitude  ratios 
for  the  vibration  of  an  unloaded  cylinder  and  for  buckling  due  to  axial  loads  are  identical.  The 

2  2 

eigenvalues  A(X,  n)  represent  the  frequency  parameter  y  oj  for  the  unloaded  cylinder  and 

2  o 

q^Q  X  for  the  buckling  problem. 

For  the  vibration  of  freely-supported  cylinders  under  axial  loads  less  than  the  buckling 

value,  the  frequency  parameter  y^  o?  and  axial  load  parameter  q9  X*  can  be  combined  to  form 

2  2  2  * 

a  single  eigenvalue  expression,  A(X,  n)  =  y  u>  +  X  .  The  eigenvalues  can  be  obtained 

without  specifying  a  value  for  axial  load.  For  this  problem  the  square  of  the  vibration  frequency 
for  any  mode  was  found  to  vary  linearly  with  axial  load. 

The  vibration  frequency  for  any  mode  and  axial  loading  can  be  obtained  by  knowing  only  the 
vibration  frequency  of  an  unloaded  cylinder,  cuo>  and  the  value  of  the  buckling  parameter  q^ 

for  that  mode.  The  relationship  between  these  parameters  is  given  by  Equation  127 

"*■  <  ( ■--£")  11271 

The  values  for  a>o  and  q2Q  can  be  obtained  from  the  figures,  where  both  values  should  be 

for  deformation  which  is  either  radial,  axial,  or  torsional.  All  vibration  frequencies  decrease 
for  a  compressive  axial  loading  and  increase  for  a  tensile  axial  loading.  For  a  given  axial 
loading  the  relative  change  in  vibration  frequency  is  greatest  for  those  modes  having  the 
lowest  buckling  values  q^. 

For  the  non-axisymmetric  modes  (n  i  0),  the  lowest  natural  vibratiou  frequencies  and 
lowest  buckling  values  for  axial  loading  for  a  given  mode  were  associated  with  deformation 
which  was  radial  except  for  long  cylinders  and  modes  with  one  circumferential  wave  where 
the  “beam”  modes  produced  the  lowest  values. 

The  lowest  vibration  frequency  for  the  unloaded  cylinder  has  one  axial  half-wave.  This  is 
generally  not  the  case  for  the  lowest  buckling  mode  for  a  cylinder  under  axial  load.  It  appears 
that  the  lowest  natural  vibration  frequency  and  the  lowest  axial  buckling  load,  for  a  given 
number  of  axial  half-waves,  have  the  same  number  of  circumferential  waves.  The  number  of 

circumferential  waves  depends  on  the  value  of  the  parameter  —■  and  the  radius  to  thickness 

Ilia 

ratio  -jj. 

The  lowest  vibration  frequencies  for  the  unloaded  cylinder  for  torsional  and  axial  vibration 
were  the  axisymmetric  modes  (n  -  0)  with  one  axial  half-wave.  The  frequencies  increase  as 
the  number  of  circumferential  waves  and  the  number  of  axial  half-waves  increase. 

The  lowest  natural  vibration  frequencies  for  the  axisymmetric  modes  (n  =  0)  are  associated 

with  vibration  which  is  radial  foi  values  of  <  2  and  pure  torsional  for  values  of  -JL  >  2. 

ma  ma 
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The  equation  for  the  uncoupled  axisymmetric  torsional  vibration  frequency  is  given  by 
Equation  121. 


2  p  (i  +v) 


m-  I,  2, 


The  only  assumption  used  to  obtain  this  equation  was  that  the  ratio  of  thickness  to  radius 
is  small  compared  to  1. 


The  lowest  axial  loading  buckling  value  for  axisymmetric  modes  also  involves  radial 

deformation  for  -f-  <  2,  and  pure  torsional  deformation  for —  >  2.  The  buckling  value  for 
ma  ma 

the  axisymmetric,  pure  torsional  mode  is  given  by  Equation  129 


P  (l-i/  ) 
q20  *  o  '  2 


(  129) 


with  the  assumption  that  the  ratio  of  thickness  to  radius  is  small  compared  to  1.  This  axial 
load,  qgQ,  is  independent  of  the  length  of  the  cylinder  and  the  number  of  axial  half-waves  in 

the  mode  shape. 


For  the  buckling  of  a  freely-supported  cylinder  due  to  a  pressure  loading,  the  axisymmetric 
modes  and  modes  with  one  circumferential  wave  have  buckling  values,  q10,  which  are  much 

higher  than  those  for  two  or  more  circumferential  waves.  For  modes  having  from  0  to  4  cir¬ 
cumferential  waves,  one  negative  buckling  value  is  obtained,  which  corresponds  to  buckling 
due  to  internal  pressure.  For  the  axisymmetric  modes  and  modes  with  one  circumferential 

Q 

wave,  a  negative  buckling  value  is  obtained  for  the  full  range  of  the  parameter-^  from  .2  to 

15.  For  modes  having  two  circumferential  wa^ss,  a  buckling  value  for  internal  pressure  exists 

only  for  values  of  <  1.  For  modes  with  three  circumferential  waves,  a  buckling  value  for 
ma  0 

internal  pressure  exists  only  for  values  of  <  .4.  For  modes  with  four  circumferential 
waves  a  buckling  value  for  internal  pressure  is  obtained  only  for  =  .2.  This  negative 
buckling  value  is  quite  large. 

For  modes  having  one  circumferential  wave,  the  lowest  buckling  value  results  for  external 

pressure  with  radial  deformation  for  values  of  — —  <  t r,  and  internal  pressure  with  a  “beam” 

ma 

mode  type  deformation  for  values  of  —  >  tr. 

J  ma 

For  the  buckling  of  a  freely-aupported  cylinder  due  to  a  constant  pressure  loading  in  an 
axisymmetric  mode  (n  =  0)  the  buckling  values  q1(j  involving  coupled  axial  and  radial  deforma¬ 
tion  are  given  by  Equation  139  for  large  values  of  . 


q  -  -  v  t  1  (139) 

10 

ft 

Numerical  calculations  for  the  range  of  radius  to  thickness  ratios  from  r-  =  20  to  5000  indi- 

a  n 

cate  that  Equation  139  can  be  used  for  values  of  >  1.  The  deformation  associated  with  these 

results  is  radial.  Both  the  positive  and  negative  buckling  loads  had  the  same  numerical  value 
for  the  ratio  of  the  radial  deformation  to  axial  deformation  and  were  of  opposite  sign. 
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For  freely-supported  cylinders  subjected  to  both  axial  loads  and  a  pressure  loading,  the 
square  of  the  vibration  frequency  for  a  given  deformation  mode  was  found  to  vary  essentially 
linearly  with  the  pressure  loading  as  well  as  the  axial  load  for  many  modes.  This  relation¬ 
ship  is  given  by  Equation  145, 


U) 


(. 


— li-) 


I20  qi0 


(14  5) 


where  wQ,  q2g»  and  q^  are  the  natural  frequency  of  the  unloaded  shell,  buckling  value  for 
axial  load,  and  buckling  value  for  a  pressure  load,  for  a  given  deformation  mode. 


The  errors  introduced  by  this  linear  relationship  are  evaluated  for  a  wide  range  of  sheli 
modes  oovering  the  full  range  of  the  shell  parameters.  Equation  145  was  found  to  give  poor 
results  for  the  axisymmetric  modes  and  modes  with  one  circumferential  wave.  For  modes 
with  two  circumferential  waves  the  results  were  within  7%  of  the  correct  values  for  all  modes 

p 

examined  for  values  of  >  .5.  The  discrepancy  decreases  as  the  number  of  eir  cumferential 

ma 

waves  in  the  mode  shape  increases.  Buckling  values  for  a  combined  axial  and  pressure  loading 
can  be  obtained  by  setting  &  equal  to  zero  in  Equation  145. 


For  those  modes  having  only  positive  buckling  values  for  external  pressure,  external 
pressure  decreases  the  vibration  frequency  and  internal  pressure  increases  the  vibration 
frequency  for  that  mode.  For  those  modes  having  negative  buckling  values  corresponding  tc 
internal  pressure,  the  vibration  frequency  appears  lo  increase  with  external  pressure  and 
decreases  with  internal  pressure  if  the  negative  buckTng  value  is  smaller  in  magnitude  than 
the  other  external  pressure  buckling  values  for  that  moce. 
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SECTION  vm 

RECOMMENDATIONS 


The  figures  contained  in  this  report  can  be  used  directly  to  obtain  the  natural  vibration 
frequencies,  buckling  values  for  axial  loading,  and  buckling  values  for  external  and  internal 
pressure  for  freely-supported  cylinders.  The  vibration  frequencies  can  be  obtained  for 
cylinders  subjected  to  axial  loads,  and  approximate  results  can  be  obtained  for  freely- 
supported  cylinders  subjected  to  both  axial  load  and  external  pre  .sure. 

Arnold  and  Warburton  (Reference  2)  suggest  the  use  of  vibration  results  for  the  freely- 
supported  cylinder  for  cylinders  with  other  boundary  conditions  by  relating  these  to  an  equiv¬ 
alent  freely-supported  cylinder.  Equivalent  wavelength  factors  were  obtained  from  experi¬ 
mental  vibration  results  for  cylinders  with  different  types  of  supports  and  are  given  in 
Reference  2.  Through  the  use  of  the  concept  of  an  equivalent  freely-supported  cylinder,  the 
figures  and  formulas  presented  coaid  be  used  to  obtain  preliminary  design  values  for  the 
vibration  frequencies  for  loaded  or  unloaded  cylinders  as  well  as  buckling  values  for  oombiaed 
loadings  for  cylinders  with  various  boundary  conditions. 

If  the  concept  of  an  equivalent  freely-supported  cylinder  is  valid  for  buckling  results  as 
well  as  for  vibration  frequencies,  the  linear  approximation  between  the  square  of  the  vibration 
frequency  and  axial  load  and  e  xternal  pressure  may  also  be  valid.  If  this  is  the  case,  vibration 
tests  could  be  conducted  in  connection  with  static  buckling  tests  to  achieve  nondestructive 
testing.  By  recording  vibration  frequencies  at  various  increments  of  load  it  may  be  possible 
to  extrapolate  to  the  buckling  values  without  damaging  costly  test  specimens  or  actual 
hardware. 

Since  the  natural  vibration  frequencies  for  cylinders  are  sensitive  to  both  boundary  condi¬ 
tions  and  imperfections,  it  may  be  possible  to  obtain  better  correlation  between  buckling  tests 
by  using  experimental  vibration  frequencies  as  an  index  of  imperfection  or  support  restraint. 
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